Photo-induced Tomonaga-Luttinger-like liquid in a Mott insulator 
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Photo-induced metallic states in a Mott insulator are studied for the half-filled, one-dimensional 
Hubbard model with the time-dependent density matrix renormalization group. An irradiation 
of strong AC fields is found to create, in the nonequilibrium steady state, a linear dispersion in 
the optical spectrum (current-current correlation) reminiscent of the Tomonaga-Luttinger liquid 
for the doped Mott insulator in equilibrium. The spin spectrum in nonequilibrium retains the des 
Cloizeaux-Pearson mode with the spin velocity differing from the charge velocity. The mechanism of 
the photocarrier-doping, along with the renormalization in the charge velocity, is analyzed in terms 
of an effective Dirac model. 

PACS numbers: 71.10.Fd, 71.30.+h, 72.20.Ht,72.40.-|-w 



Introduction — Doped one-dimensional(ld) Mott in- 
sulators are fascinating due to a special metallic state 
known as the Tomonaga-Luttinger (TL) liquid, where ex- 
citations have collective nature as distinct from the con- 
ventional Fermi liquid 0- Specifically, the charge velocity 
becomes renormalized due to the electron-electron inter- 
action, so that the charge excitation propagates with a 
velocity different from that of spin excitations — a hall- 
mark of the "spin-charge separation" . Experimentally, 
TL liquids have been observed in quantum wires and car- 
bon nanotubesQ. 

Now, there is another way of making a Mott insula- 
tor metallic, which is entirely different from the chemical 
doping. Namely, photo-doping is now being highlighted 
as a way to control the carrier density in Id correlated 
systems II, H, H S 11 H [iS- Pump-probe experiments 
have shown that strong electric fields can turn Mott in- 
sulating crystals into metals [1, 0, 0]. Theoretically, the 
change of the conductivity by irradiation was studied in 
with exact diagonalization. Features in the photo- 
doping conceptually distinct from the chemical doping 
are: (a) the system is out of equilibrium, and (b) two 
types of carriers, i.e., electrons (=doubly occupied sites) 
and holes, coexist because they are pair-produced. Al- 
though electron-hole systems in equilibrium have been 
studied in the past in terms of the TL theory llj , prop- 
erties, especially the spectral properties, of the nonequi- 
librium metallic states induced by irradiation are yet to 
be understood. 

This has motivated the present work, where we shed 
light to this problem by studying the optical (current- 
current) and spin correlation functions for a nonequilib- 
rium steady state in Id Mott insulators in strong AC elec- 
tric fields. This is done in two steps: we first employ nu- 
merical simulations with the time-dependent density ma- 
trix renormalization group (td-DMRG) to show that 
we have a TL-like linear dispersion in the photo-doped 
system, where the charge and spin excitations have dif- 
ferent velocities. We then confirm the result from a field 
theoretic result. A starting point is ref.fH, where 
the present authors pointed out that there is an inter- 
esting analogy between the Schwinger mechanism in the 



decay of the quantum electrodynamics (QED) vacuum 
governed by the electron-positron creation ratefisj. and 
the dielectric breakdown of Mott insulators in DC elec- 
tric fields governed by the electron-hole creation rate. In 
Id, correlated electron systems can be mapped to a sim- 
ple Dirac model via the massive Thirring model, which 
is dual to the sine-Gordon model for charge degrees of 
freedom of the Id Hubbard model. We then analyze the 
photo-doping with the effective Dirac model. There, the 
optical spectrum is formulated with the Floquet analy- 
sis for strong AC external fields. The result reproduces 
the main features in the td-DMRG result, providing a 
physical picture for the dynamics of the nonequilibrium 
collective modes in the photo-induced state. 

Nonequilibrium steady state — We consider a Mott 
insulator in strong AC electric fields in the half-filled. 
Id Hubbard model with the Hamiltonian given, in stan- 
dard notation, by H{t) = Hq + Hp(t), where Hq = 



^hop X]i(T ('^l+lo-'^^o- 

Fe{t) sin(m) 



'it 



and i?F(i) = 
rill). Here 



F and Q are, respectively, the strength and frequency of 
the external electric field, which is switched on at t = 
(hence the insertion of a step function). In the calcula- 
tion we take the length of the system L = 80, a time step 
At = 0.04, and the DMRG Hilbert space size of m = 140, 
in natural units. After obtaining the groundstate |^'o) 
of H{t < 0) with the finite-system DMRG algorithm 
p^ . we let the system evolve according to the time- 
dependent Hamiltonian H{t) with the td-DMRG [H] to 
obtain the wave function |^'(<)) = C/(t; 0)|^'o), with the 
time-evolution operator U{t;t') = fe~'^t' (f : the 

time-ordering). 

For < > 0, the system relaxes into a 
nonequilibrium steady state, where we can 
define the photo-doping rate as Xph{t) = 

rEi(avg(*(i)|n»T"al*W) " (^'o|'^»T"a l*o)) which 
is the increment in the double occupancy. We also 
monitor the total energy, Ea.-vg{t) = avg(*(t)|i?(t)|*(t)). 
In these expressions, we eliminate the Tperiod = '2,1^/^ 
oscillation by taking the average over each period, as 
denoted by "avg" . The time profiles of the photo-doping 
rate and the total energy in FigU] are similar, which 
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FIG. 1: (Color online) Inset: Time-dependent DMRG result 
for the total energy(a), photo-doping rate(b) for the half-filled 
Id Hubbard model with C//thop = 8.0, ^ = 8.0. The main 
panels depict their final values at t = Ti vs the field strength 
F, where the solid lines are parabola fitted to the data. 

indicates that the energy absorbed from the AC field 
is used to excite electron- hole pairs, i.e., photo-doping. 
The system relaxes to a steady state until the situation 
where pair production rate = annihilation (stimulated 
emission) rate is achieved. 

The doping rate depends both on the frequency and 
the strength F of the electric field. When Vl is greater 
than the Mott gap A, excitations occur via one-photon 
absorptions. Even when f2 < A, multi-photon processes 
can excite the system above the gap for sufficiently strong 
fields, where the process becomes resonant when mfi = A 
with m being an integer [I^ . The extreme case is the DC 
limit f2 — > 0, where many-body Landau-Zener tunneling 
across the Mott gap induces metallization [H, [3, [13 • 

Correlation functions — To characterize the Id many- 
body system we calculate the correlation functions after 
a steady state is attained, 




FIG. 2: (Color online) Current-current correlation 
Imxjj(f,Ti;i, j)(a), the local spin energy eP'"(t)(b), and the 
double occupancy nf(t)(c) for zero field _F = (left pan- 
els), and for a finite field F = 0.2 (right) with U/thop = 
8.0, r^/fhop = 8.0. Here the probed position j in the pump- 
probe process is set to be at the center {j = 40) of the fd 
system of length 80. Red broken lines indicate the region for 
relaxation due to spin (see text). The small ripples (c, right 
panel) is the T = 2n/Q period fluctuation. 



XAB{t,Ti;i,j) = {'f{t)\A,U{t,Ti)B,\'f{Ti)), (1) 

where A, B are operators, e.g., the current 
Ji = -ithopJ2cri4+ia(^i'y - h.c). Or the spin 
Si = ^J2ap4a^ai3Cii3, and Ti is the time around 
which the steady state is reached (and the curves in 
FiglD fiatten; typically Ti = 50, 100 depending on the 
field strength). Physically, XJJ represents the probing 
process in standard pump-probe experiments, where a 
photon in the probe light generates a local electron-hole 
pair at position j. To focus on the interplay between 
the charge and spin degrees of freedom, we can compare 
in Fig. [2] the current correlation function with the 
behaviors of spin and charge. For the spin we define the 
local spin energy, er"(t)/ Jexch = (^'•^(i)|s,+i • 
which measures the exchange energy (normalized 
by the exchange coupHng Jcxch = 4f^pp /£/)f20j. 
For the charge we examine the double occupancy 
defined by nf{t) = («'"'(t)|ni|n,x|«'"'(<)). Here 
\-^''{t)) ^ C/(<,Ti)Jj|*(Ti)) is the state where a 
perturbation (Jj) is added at site j on the steady, 
nonequilibrium state. The behavior of the three 
quantities shows that the temporal evolution after the 
probe-excitation on j a,t t = Ti propagates in two 



processes. The first is diffusion of the doublon-hole 
pair, which is followed by the relaxation process. The 
relaxation is seen as a decay of the current correlation, 
which is seen to be accompanied by a disturbance in the 
spin structure (as marked with a red broken line in the 
figure). This indicates that the spin and charge degrees 
of freedom become coupled more strongly, where spins 
act as a kind of energy reservoir for charges. Spin-charge 
coupling already exists in equilibrium for higher-energy 



states. |2l| which is natural since charge excitations act 
as boundary conditions to spins with spin vanishing 
at doubly occupied or empty sites. The decay of the 
current correlation, already present for zero AC electric 
field {F = 0), becomes faster in finite fields, which 
implies that the spin-charge coupling becomes stronger, 
due to higher-energy states become involved. However, 
in the field range studied here, the coupling is not strong 
enough to destroy the spin-charge separation picture 
[l0| . i.e., the spin and charge degrees of freedoms still 
have independent dispersions, as discussed below. 

Collective excitations in nonequilibrium — So what is 
the nature of the photo-induced carriers? We can obtain 
the excitation spectrum as the Fourier transform of the 
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FIG. 3: (Color online) The optical spectrum, lmxjj{q,u!), 
with ~ 8.0 for (a) the half-filled case with no AC fields, 

(b) the half-filled with a finite AC field F/thop = 0.1, and (c) 
the electron-doped with no AC fields {n = 1.025, F — 0). 



correlation functions, 



(2) 



,iu]{t-Ti)-iq(j-jc) 



lmXAB{t,Ti;j,jc 



where T2 is set to be prior to the time at which the wave- 
front reaches the sample boundary. We call the quantity, 
following the equilibrium case, the optical spectrum for 
Xjj, and the spin spectrum for Xss- 

A striking feature in the result in Fig. [31 a first key 
finding here, is that, while we have an optical gap (= 
Mott gap) in the zero AC field (Fig. [3)Ja)), a set of new 
states with a zero-gap, linear dispersion, ujq ^ i'ptl|'?|, 
emerges in the gap in finite AC fields (Fig. [3]^b),(c)). 
We call the metallic state a photo-induced Tomonaga- 
Luttinger-like liquid in the sense that its charge veloc- 
ity wpTL is renormalized by the electron interaction. Its 
value is in fact similar to the equilibrium counterpart as 
can be seen by comparing the slope in nonequilibrium 
with that for the equilibrium, doped system (Fig. EJ^e)). 

If we turn to the spin spectrum in Fig. 31 we can 
clearly see the des Cloizeaux-Pearson mode in equilib- 
rium (precisely speaking, this was obtained by the Bethe 
ansatz method |23|). which survives in finite fields . We 
note that the charge and spin velocities defined and nu- 
merically obtained in nonequilibrium here are different. 
So in this particular sense we have a spin-charge sepa- 
ration. As the field becomes stronger, the spectrum be- 
comes blurred. Specifically, the antiferromagnetic fluctu- 
ation represented by the peak value Xssiq 7r,(^ ~ 0), 
becomes smaller, which should be because the induced 
carriers act to melt the magnetic order. 

An effective model — Due to its many-body nature, 
it is difficult to obtain a mathematically rigorous anal- 
ysis of the photo-induced metallic state in Mott insu- 
lators. However, we can carve out some of its proper- 
ties in terms of an effective model, which was initiated 




FIG. 4: (Color online) (a) Spin correlation lmxss{t,Ti;i,j) 
for zero field (left) and for a finite field F = 0.4 (right) 
with U/thop = 8.0, Q/thop = 8.0. (b) Color-coded spin 
spectrum ImXss(g,'i') in the half-filled Hubbard model with 
= 8.0 for zero AC field(left), or for a finite AC field 
F/t\iap = 0.4(right). The white dashed lines are the des 
Cloiseaux-Pearson mode, (c) The AF peak \mxss{ci,uj = 0) 
vs q are also displayed for various field strengths. 



by Luther, Emery, and by Giamarchi to study linear- 
responses in the presence of a charge gap [1^ Hj, [2^, . 
In this approach we start from the 1+1 dimensional mas- 
sive Thirring model. The charge degrees of freedom is 
represented by two spinless fermions ^' = {ipi{x), ^2{x)), 
where ^pi,^p2 represent left mover and right mover, re- 
spectively, and the Hamiltonian reads 



Hf'{t)^v, dx 



*T(x)(-iV,(0a3 + -ai)*(x) 



Hii?,) 



where <Ti is Pauli matrices, Vc the renormalized charge ve- 
locity, and A the Umklapp-scattering coupling constant, 
an ascendant of the original Mott gap at half filling. This 
model is dual to the sine-Gordon model, where the size 
of the interaction term Hj = g J dx[{^^^)^ - (^ViVl')^ 
translates to the sine-Gordon coupling through a duality 
relation[23|. However, in the following we make a further 
simplification, namely, we neglect the self-interaction 
Hj. This by no means implies a neglect of the origi- 
nal electron-electron interaction, but amounts to neglect 
self-energy corrections to quasi-particlc life time, etc. We 
employ this approximation to focus on the effect of the 
charge gap, which primarily appears from the first term 
in eqn.Q. 

The AC electric field is taken as the coupling, V^, {t) = 
dx + iAi{t), in the above, with Ai{t) = {F/fl)smQt. 
We then obtain the nonlinear, nonequilibrium evolu- 
tion of the state in this model. After a Fourier trans- 
form, the equation of motion becomes i^|4'fc(t)) ~ 
[vc{k + (F/il) sin nt)a3 + y'^i)] Here we adopt 

the Floquet method for treating systems in AC fields 
(see e.g., [l^), i.e., we seek a solution of the form 



-iSa {k)t—imQ,t 



u^{k)), where m is the 
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FIG. 5: (Color online) (a) The optical spectrum 
Imx''"^'"^('i', of the effective Dirac model in a finite AC field 
with f2/A = 1.5, F = 0.1. (b) The Floquet quasi-energy 
Eaiq), where the red lines represent occupied states, green 
lines unoccupied states, (c) A schematic level repulsion. 



Floquet index, a labels the eigenstate, and e„ is Flo- 
quet's quasi-energy. From the equation of motion, 
the Floquet modes satisfy a set of linear relations 
En(^eff)™"K(fc)) - ea(fc)K(fc)), where the effective 
Floquet Hamiltonian has matrix elements (i/eff)™™ = 



m.mibl 



TiVcjyCTs. We con- 



sider a sudden switch on of the AC field, i.e., the time evo- 
lution starts from the groundstate |V'o(^)) of Vckas+Y^^i- 
The solution is then \^Q{k;t)) = '^^(f>a{k)\ua{k;t)) with 
4>a{k) = (u^(k)\tbo(k)). The current correlation function 
in this model is [29[ 

Imx°'-^(9,-) = 2e^| g 

X (1^/3^=)^ ~ \Mk + q)\^)S{ep{k) - e„(fc + q) + uj) 

xTr[K(fc + q))«"(fc + g)|r3|u^(fc))(u^(fc)|r3] . (4) 

The optical spectrum, FiglSlJa), shows that a metallic, 
linear-dispersion mode emerges in the gap in the pres- 
ence of an AC field. The result does resemble the present 
td-DMRG result (FiglS]). The spectrum, despite being 
nonequilibrium, is also similar to the massive Thirring 



model result for the doped Id Mott insulator. |23j The 
origin of the gapless excitation in nonequilibrium can be 
traced back, in the present massive Thirring 4- Floquet 
analysis, to the quasi-energy level scheme in Fig. [5l^b). 
In the absence of AC fields, the Dirac model has two 
branches with the lower band completely filled. As we 
turn on the AC field, we have a series of replicas (i.e., Flo- 
quet modes equally spaced by fl) for each of the electron 
and hole branches. Physically, these modes correspond 
to m-photon absorbed states. As in standard quantum 
mechanics, level repulsion takes place when two modes 
cross with each other with nonzero matrix element be- 
tween them. In our model, the most important feature 
is the level repulsion between an occupied Floquet mode 
(red in FiglSI^b)) and an unoccupied one (green). As 
shown in the blowup (Figl^Kc)), gapless excitations then 
emerge across the occupied and unoccupied states, which 
should contribute to the optical spectrum when we eval- 
uate eqn.([l]). Specifically, the linear dispersion of the 
photo-induced state is ~ Vc\q\, where the renormalized 
velocity v^, a parameter independent of the spin velocity 
Vs cx t'^ap/U, is conceived to have a value that depends 
on the detail of the irradiation. 

Another interesting point is when the photon energy 
is below the gap A. In this situation, multi-photon 
processes with m > 2 are necessary to excite the system. 
We can indeed show that, when the condition U ^ A/m 
is met, the excitation becomes resonant and carriers are 
injected efficiently, which will be described elsewhere. 

In conclusion, we have shown that a Mott insulator 
irradiated by strong AC electric fields has a collective 
mode which is reminiscent of the Tomonaga-Luttinger 
liquid in equilibrium. While spins and charges are cou- 
pled in the relaxation process, the collective modes have 
different spin and charge velocities. Emergence of the lin- 
ear collective mode is also supported by an effective Dirac 
model. We have greatly benefited from discussions with 
Naoto Tsuji on the Floquet method. TO acknowledges 
Keiji Saito and Keiichiro Nasu for valuable discussions. 
This work has been supported in part by a Grant-in- Aid 
for Scientific Research on a Priority Area "Anomalous 
quantum materials" from the Japanese Ministry of Edu- 
cation. 
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